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ABSTRACT 

We consider the electrodynamics of in-spiraling binary pulsars, showing that there are two 
distinct ways in which they may emit radiation. On the one hand, even if the pulsars do 
not rotate, we show that in vacuo orbital rotation generates magnetic quadrupole emission, 
which, in the late stages of the binary evolution becomes nearly as effective as magnetic 
dipole emission by a millisecond pulsar. On the other hand, we show that interactions of the 
two magnetic helds generate powerful induction electric helds, which cannot be screened by 
a suitable distribution of charges and currents like they are in isolated pulsars. We compute 
approximate electromotive forces for this case. 
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1 INTRODUCTION 

The recent discovery of the binary pulsar J0737-3039A/B jBurgay| 
|et al.|2003] |Lyne et al.|2004^ has highlighted the fact that pulsars 
may retain their magnetic fields for at least a significant fraction of 
the time it takes them to in-spiral to their eventual merger. 

This discovery has been used so far to illustrate the extraordi¬ 
nary opportunity to study hitherto unaccessible General Relativis¬ 
tic effects, while less attention has been paid to the electromagnetic 
interactions of the two magnetospheres (but see |Kramer & StaTrs| 
|2008| also for a discussion of magnetospheric effects). 

But even before this discovery, |Vietri|fl996| > and |Hansen & L yu^ 
|tikov| ( [200r| l have discussed some of the consequences of assuming 
that the two inspiralling, merging pulsars, (i.e., not just neutron 
stars) may produce observationally interesting electromagnetic sig¬ 
nals. In particular, [Hansen & Lyutikov| ( |200l1 considered the mag¬ 
netospheric interaction of two inspiralling pulsars with very differ¬ 
ent magnetic fields: a fast, but weak-field pulsar, and a slow mag- 
netar. 

In this paper, we wish to take a different point of view, that in 
which the two fields do not differ by four to six orders of magnitude 
like in [Hansen & Lyutikovj pOOT) , but are more evenly matched, 
though not necessarily of the same order of magnitude. In particu¬ 
lar, we wish to stress first that the ultimate energy reservoir to be 
tapped is the binary orbital motion, which we do by considering 
the simplest model possible (two point-like dipoles orbiting their 
common center of mass in vacuo), and pointing ou that this sim¬ 
ple estimate alone leads to an observationally interesting signal. 
The reason is that binary pulsars will rotate ever faster as a con¬ 
sequence of the orbital decay induced by gravitational radiation. If 
the pulsars manage to retain some signifiant fraction of their mag- 

* email: emanuele.sobacchi@sns.it 
f email: mario.vietri@sns.it 


netic fields until the moment of merging, their evolution will lead 
naturally to a sub-millisecond object emitting copious amounts of 
radiation in the last instants before merging, even when we neglect 
the pulsars’ rotation around their spin axes. We shall first discuss 
this emission mechanism, which occurs in vacuo. 

On the other hand, under these circumstances there is a natural 
mechanism that will create an induction electric field with a compo¬ 
nent directed along the magnetic field, for an arbitrary orientation 
of the axes in question. This mechanism follows from the well- 
known inability of the two pulsars to synchronize their spin periods 
with their orbital period, a direct consequence of their tiny sizes (in 
terms of volume, not mass, of course; see |Bildsten & Cutler[1992| >. 
Since sufficiently late in the binary life the orbital period is much 
shorter than the spin periods, we may picture the two pulsars as 
non-rotating; still, the orbital motion swings each magnetosphere 
by the other one, so that any particular part of it is compressed on 
the day-side, and is free to expand on the night-side. This magne¬ 
tospheric pumping implies a time variability of the local magnetic 
field, and in turn this leads to an induction electric field. 

This situation differs considerably from the aligned rotator (|Con-| 
[topoulos et al.[[T999[ l, where a judicious choice of a stationary 
charge distribution inside the magnetosphere can short out the com¬ 
ponent of the electric field along the magnetic field; we shall show 
later that this mechanism cannot be effectively screened, as it is in 
the case of an isolated, rotating aligned pulsar (IRAP for short). 
This suggests that large-scale, large-amplitude electric fields not 
necessarily orthogonal to the magnetic field do exist in the magne¬ 
tospheres of binary pulsars. 

It is the aim of this paper to investigate some consequences of 
this simple idea for the (photonic) observability of in-spiraling bi¬ 
nary pulsars as GWR pushes them closer and closer. In the next 
Section, we will derive the energy lost per unit time by two co¬ 
rotating magnetic moments in vacuo, a generalization of [Pacini] 
l[1967[’s formula for the isolated magnetic moment. In Section 3, 
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we consider in some detail the induction effect briefly discussed 
above, still using the in vacuo case as our model, for lack of a fully 
satisfactory realistic model. We shall discuss how the component 
of the induction electric field along a magnetic field line cannot be 
shorted by arbitrary charges and currents located in the magneto¬ 
sphere, and we shall provide some exact and some approximate es¬ 
timate of the available electromotive forces. The last Section sum¬ 
marizes our results. 


2 IN VACUO 

The starting point of this paper is a computation of the amount of 
energy radiated per unit time by a configuration of two point-like 
dipoles fli, p ,2 rotating around their common center of mass. We 
assume in this that the stars carrying the magnetic moments are 
not spinning at all (fl = 0 ), and call cu 7 ^ 0 the orbital angular 
frequency. 

This configuration does not emit magnetic dipole radiation, be¬ 
cause the total magnetic moment / 2 i + jl 2 is a constant, thus to this 
order an observer located at infinity will perceive no time-varying 
field. However, said observer will perceive a time-varying field to 
higher order, due to the fact that the magnetic field is stronger when 
the star with the stronger dipole moment is closer, and weaker half 
an orbital period later when it is furthest, and the weaker dipole is 
closest. 

This effect depends on the orbital radius r: the modulation of the 
fields vanishes for r = 0 , and it is clearly linear in r for r <C the 
observer’s distance. Thus the radiated power will be quadratic in r, 
showing that we are dealing with quadrupole magnetic radiation. 

It is obviously possible to obtain the radiated power by means 
of a vector-spherical harmonics analysis ( |Jackson|1975| l. However 
it is also possible to reach the same result via a simpler approach, 
which also allows us to introduce the fields we shall use in the next 
Section. 

|Monagh^ ( |1968^ has given exact expressions for the electro¬ 
magnetic potentials and fields generated by an arbitrarily moving 
dipole, i.e., a particle endowed with both electric (jp) and magnetic 
(jl) dipoles. We just need to specialize to the case when, in the 
particle (actually a star, in our case) reference frame, it has only 
a non-vanishing magnetic dipole moment; in this case Monaghan 
shows that, when the pure dipole magnetic moment is moving with 
arbitrary speed c/3, it appears to have an electric dipole electric mo¬ 
ment 


p = /3 A // . 


( 1 ) 


In this case the electromagnetic potentials are: 

X _ ( P ^ , rf/rAn ^ d P A fl\ 

^ [ KR^ ^ ~dt KRc ^ ~dt KRc I 

\ / ret 


f A p) ■ n d 0Afl)-n \ 

\ KR^ dt KRc I 

\ / ret 


where K — 1 — n-f5,fi = v/c and jj is the magnetic dipole of the 
pulsar. 

|Monagh^ (| 1 968|l also gives an especially useful expression for 


the actual fields: 

3(p • h)h — p 
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c dt 
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where of course eq.[T]applies. Correspondingly, one also obtains: 
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The reason why these expressions are useful is that all radiation 
terms are contained in the last term on the rhs of eq.|^(an d 0 : 

j 2 


Ersid - r. 


1 d^ fhA{hAp) — pAh 


dt^ 


KR 


( 6 ) 


Like we said above, we completely neglect the star’s rotation, so 
that pis a constant and p = PAp. Expanding the above expression 
to lowest order in /3 we find: 


p _ {n- p.)0Ah) 


(7) 


We can now specialize to the case of a binary pulsar, each with 
magnetic moments pi and p 2 ,hi = n 2 , and /3i = —p 2 = P (for 
equal mass stars, of course) to obtain 


Erad - 


(n- {pi -p 2 ))(p Ah) 
Rc^ 


Br3,d - Tl A Ej-: 


( 8 ) 


From Poyinting’s vector, cEmd A Bm^d 14:71 we obtain the total ra¬ 
diated power as 


f ^ {n-iPi- P^nPAnf 
J 47rc'^ 

We can now, first average over the orbital phase, then carry out the 
integral over angles (see the Appendix) to obtain: 

^4/q2 

Pm ={3{iii - P2)\ + 4{iii - pL2)l) ( 10 ) 

where pz the component of the magnetic field along the orbit’s nor¬ 
mal, while pL± is the component in the plane of the orbit. 

This formula clearly describes the radiation from a magnetic 
quadrupole, it is a full factor P"^ smaller than the dipole term. It has 
partially been derived before: [Harrison & Tademaru| ( |1975) stud¬ 
ied the radiation from a pulsar with a (single) magnetic dipole dis¬ 
placed by a fixed amount from the rotation axis. In their case, they 
had dipolar radiation because the components of p in the orbital 
plane do not remain constant as a consequence of the star rotation. 
However, the component produces a quadrupole term which 
they computed, and coincides with our result above. 

We now specialize to the case of a binary pulsar whose orbital 
decay is driven by gravitational wave radiation losses, in which case 
(as we are about to see) losses are completely dominated by GW. 
It can be seen from the equation above that electromagnetic losses 
depend on the star-to-star distance a as a“^, and are thus strongly 
peaked around the first moment of physical contact; we may thus 
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safely assume that the orbit has already been circularized by GWs, 
and consider only circular orbits. We have 


2 


OJ 


GM 


E = 


1 GMM 

2 a 


( 11 ) 


where the total mass and reduced mass are given, respectively, by 
M = Ml + M 2 , M = M 1 M 2 /M; the GW loss rate is 


Egw 


From the above we find 


da 

dt 


32 G* M^M^ 
5 c® a® 

2 Egw a? 
GMM ' 


( 12 ) 


(13) 


Thanks to the equation above, in the limit Mi = M 2 , we can trans¬ 
form eq.|10|into 

dEyi ^ 1 jdx 
da 384 


where we have called Em the amount of orbital energy lost via 
electro-magnetic processes (Pm = dEm / dt) and 

_ 3(/ii —/22)i + 4(/riz — p2z)^ 

x= - ——2 - (15) 

MI +K 


and fd = Hi + ^ 2 - Integrating over a we obtain the total amount 
of radiation emitted by the magnetic quadrupole term: 


Pm 


1 




1152 at 


(16) 


For typical values fi = 10^° G cm®, x « 3 and Omin = 2 x 10® cm, 
we obtain 


Pm 


6.5 X 10® 


M 


1030 G cm® 


2 

erg . 


(17) 


This amount of radiation must accompany the final stages of 
a binary pulsar. It is small compared to the total amount radiated 
away by GWs, as we anticipated. It is strongly concentrated to¬ 
ward the latest moments in the binary existence, because Pm oc 
oc The above estimate assumes that fields as large 

10®® G exist up to the time of merger. However, even a field as low 
as 10®® G may produce a detectable amount of radiation. Thus the 
above mechanism essentially predicts the existence of radio bursts 
of roughly millisecond duration. 

Though admittedly the estimate above is a text-book exercise, 
we have been unable to find it in the literature. Its importance lies 
in the illustration of the fact that the ultimate energy reservoir is 
the binary’s orbital energy, and in the fact that, even under highly 
idealized assumptions, it leads to a potentially detectable signal. 

Still, we show below that another charge-acceleration mecha¬ 
nism exists, which is likely to lead to an estimate of radiated elec¬ 
tromagnetic energy dwarfing this one. 


3 INDUCTION 

In order to see why induction electric fields will be present when 
both stars have a sizable magnetic field, we consider first two stars 
in vacuo, neglecting the stars’ rotation but including their orbital 
motion. In this case, the side of the field which is brought by or¬ 
bital motion closer to the companion is compressed by the compan¬ 
ion’s magnetic field, but, as theorbital motion leads it away from the 
companion, it is free to expand again. It is this rhythmic compres¬ 
sion which generates locally a transient magnetic field, hence an 
induction electric field. 


This periodic effect is illustrated in Fig.[2 where the position of a 
magnetic line has been computed for identical, point-like magnetic 
dipoles. 

This effect exists also in the opposite limit, when the orbital pe¬ 
riod is much longer than at least one of the spin periods: again, 
when a given portion of the magnetosphere is on the day-side it is 
compressed by the companion’s field, while it expands as it moves 
toward the night-side. 

There are two special circumstances under which this effect does 
not take place. If the stars could synchronize their orbital and rota¬ 
tional motions, thanks for instance to tidal forces, then they would 
present to their companion always the same portion of their mag¬ 
netic field, which would not be rhythmically compressed. However, 
following the work of |Bildsten & Cutler|jl992^ , we know that neu¬ 
tron stars are too small for tidal forces to produce orbit-spin syn¬ 
chronization. 

Also, the induction field vanishes if the two spin axes and the 
orbital angular momentum are parallel, as a consequence of the 
system’s reflection symmetry. 

We now consider what happens when we allow for the presence 
of free charges inside the magnetosphere. At this point, one might 
reason as |Goldreich & Julian| ( |196^ : given the abundance of free 
charges, and the fact that they move freely along magnetic field 
lines, they may redistribute themselves so as to short out the com¬ 
ponent of the electric field along B, yielding the transversality con¬ 
dition. 


E-B = 0. 


(18) 


This hunch is fully home out by the numerical solution of |Con-| 
|topoulos et al.| ( [T99^ in the (stationary!) case of the isolated, rotat¬ 
ing aligned pulsar (IRAP for short). In a time-dependent situation, 
the same spontaneous redistribution of charges will try once again 
to enforce the transversality condition, eq.[T^ and it will succeed to 
the extent that we may assume charge carriers to be infinitely fast, 
but we may safely expect this to fail, to some extent, when account 
is taken of the finite speed of charge carriers. 

This situation is similar to the penetration of electromagnetic 
waves into a good conductor, a textbook problem ( |Jackson|1975[ ): 
so long as the speed of the charge carriers is assumed infinite, there 
can be no penetration of the wave fields inside the conductor, but 
when account is taken of either the finite speed of the charge carri¬ 
ers, or of the existence of a non-zero, albeit small, resistivity, then 
the wave fields are able to penetrate for a skin depth inside the con¬ 
ductor. 

In our problem, we do not really have to take account of a change 
of speed for charge carriers due to a finite resistivity: after all, mag- 
netospheric electrons and positrons are likely to be always fully 
relativistic, hence 1 — u/c <C 1. Instead, we must take into account 
the fact that, as GWs drive the binary into ever closer orbits, the 
forcing-perturbation time-scale becomes comparable to the time 
it takes perturbations to propagate inside the magnetosphere, and 
thus to the time it takes to restore equipotentiality. This way the 
assumption that charge carriers are moving extremely fast is surely 
violated. 


We may derive the orbital time from eq. 11 using r = a/2 


as each star’s distance from the center of mass of the system, 
Torh = 4v/27rr®^®/(GM)®'^®.Forperturbationspropagatinginthe 
poloidal plane, we may take as the time to re-establish equilibrium 
the time T^c that charge carriers take to follow (at speed c) a pure 
magnetic dipole field line; these lines have equation d = D sin® 6, 
where the constant D is the maximum distance of the line from the 
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Figure 1. Deformation of field lines from dayside (red) to nightside (green) for two identical, point-like dipoles in (—r; 0) and (r; 0) oriented (respectively) 
at 90°and 45°; the line in the upper panel, if unperturbed, is closer to the surface of the star. 


star, occurring on the equatorial plane ^ = 7r/2. Its length is 


L= [\{6)d6='^D , (19) 

io ^ 

hence Tcc{D) = 7r_D/(2c). For the mid-distance line, D = r, we 
find: 


X 


= r) ^ J_ /GM\ 
Torh 8%/2 V rc'^ ) 


0.057 


/ M 10 km i?NS 
\2.8 X Mq i?NS r 


1/2 


( 20 ) 


For perturbations propagating along the toroidal direction, we can 
directly compare the perturbation speed (which equals the orbital 
speed Vorb) to the speed of light, finding 


/ 

X 


‘t’orb 

C 


0.22 

1 ^^2 ; 


0.22 


M 10 km i7NS \ 
2.8 X Mq 77ns J 


1/2 


( 21 ) 


From these equations, we see that at large distances x, 'X ^ 
hence we expect an accurate screening to occur, but as the two stars 
approach, x ~ 0.1, x^ ~ 0.2, which means the screening will be 
less than perfect. As a comparison, in a metallic onductor cr rs 
5 X 10^^ while, for optical wavelengths, u ~ 5 x 10^^ Hz, 
giving X = ~ 10“^. 

By analogy with the finite-resistivity conductor mentioned 
above, this implies that perturbation fields, due to the time- 
dependent perturbation caused by the companion star, will pene¬ 
trate for a finite length inside the (otherwise shielded) magneto¬ 
sphere, without the local charge carriers being able to short out the 
component of E along B, hence 

E- B^O . (22) 


We can actually push the analogy further to discuss the skin depth. 
In conducting media with conductivity a (in planar geometry!), the 


wavenumber is found to be ( |Jackson|1975t : 

kx2Tv(l + t)y^ (23) 

c 

where v is the impinging wave frequency. In our case, we can take 
= 1/Torb and a « l/Tcc, because I/ct is the time scale over 
which a local charge excess spreads itself over the surface of a con¬ 
ductor, i.e., the time scale to restore equipotentiality, which is the 
same physical interpretation of Tec- We obtain for the skin depth 

S = c{TeeTarby^^ ■ (24) 


As a sanity check, we check that for very large separations (r —>■ 
oo) shielding is found to be effective, as expected. In this limit, 
we cannot neglect the stellar rotation rate with respect to the 
orbital one, lj, because Q ^ ljj; each star will then have its own 
corotating magnetosphere out to a distance c/fl, with typical time 
Tcc(c/H) = 7r/2f2 (or shorter for inner lines); in this case, Tcc does 
not depend on r. Tab oc and the ratio S/r oc —>■ 0, 

showing that (at large separations) the skin depth is exactly that, a 
thin layer where perturbations can penetrate. It is worth remarking 
that, in the case of the binary pulsar, (5/r « 4.8 for rcc(c/fl) for 
PSR J0737-3039A, and a fortiori 5/r > 1 for the slower PSR 
J0737-3039B. 

Conversely, when r —> i?NS, 


S 

r 


cTeejr) - 1/2 _ TT 


2x1/2 


= 0 ( 1 ) 


(25) 


which shows that the skin depth is as large as the whole magne¬ 
tosphere. In this limit, the analysis leading to the expressions for 
k and 5 becomes invalid. In fact, objects of spherical symmetry 
smaller than about the mean free path for deflection of conduction 
electrons (i.e., in our case, cTcc) are generally considered to be only 
weakly affected by skin depth effects, in the sense that the spherical 
metal conductor is not shielded from outside fields, see for instance 
|Petrov| ( |l981| l; |Morokhov et al.| ( |1981t . For this reason, and because 
5 « r at the moment of merger, we feel justified in neglecting skin 
depth effects in the latest stages of the binary evolution. 
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Figure 2. The total electromotive force along a magnetic field line closing at 
the midpoint between two stars of equal magnetic moments, as the distance 
between the stars is reduced by GWR. The solid line is an fit. 

Given the arbitrary orientation of the three axes involved, 
I, fl\,pL 2 , and the fact that the skin depth, in the late stages of the bi¬ 
nary evolution, is of order of the orbital separation, <5 « r, it seems 
difficult to escape the conclusion that, under these circumstances, 

S • B / 0 . (26) 

3.1 Numerical estimates 

We expect emission by particles accelerated by (unscreened) in¬ 
duction electric fields to be strongly concentrated toward the last 
moments before the binary merges, exactly like the quadrupole ra¬ 
diation we considered above. Hence we shall neglect shielding, as 
justified above. Also, we are interested in the total electromotive 
force along closed magnetic field lines, both because we still ex¬ 
pect synchrotron losses to constrain particles to follow magnetic 
field lines, and because the corotating magnetosphere is where the 
largest particle density is located, and is thus likely to be where 
most radiation is emitted. 

In order to get an idea of the size of this effect, we use as a 
model the in vacuo magnetosphere, using once again Monhagan’s 
expressions, since they represent the exact solutions for the fields 
of point-like magnetic dipole, in the in vacuo case. We choose arbi¬ 
trary orientations of the magnetic dipoles such that (/2i A/ 22 ) • ( 7 ^ 0, 
where I is the orbital angular momentum. A configuration for the 
magnetic field for a given, generic choice for the magnetic moments 
is given in fig.[T] 

On the line connecting the two dipoles, the two magnetospheres 
touch approximately where ftt/Ri ~ with Ri + R 2 the 

stellar separation. We now compute numerically, for the in vacuo 
case, the integral 

£ = J E-ds (27) 

for the part (outside the star) of a magnetic field line which passes 
close to Ri, i.e., roughly the line separating the two magneto¬ 
spheres. The numerical result for the exact fields given by Mon- 
hagan (eqs. |^and|^ is displayed in fig.|^ Here we have chosen 
fii = /r 2 = 10^° in cgs units, from which i?i « R 2 (the approx¬ 
imate sign, instead of the equality sign, results from the different 
inclinations of two dipoles). Please notice that the different dots 
refer to different field lines: we are comparing the electromotive 
force always for the line separating the two magnetospheres, which 
changes as GW reduces the stars’ separation 2r. 


Figure 3. The electromotive force for lines reaching only as far as d, at a 
given instant. The solid line is a approximation. 

Our approximation, £ oc is also displayed in fig.|^ It can 

be obtained by means of the following argument. We are interested 
in £ when the stars are very close and in a region very close to 
both stars, in which case we see from eqs.|^and[^that the fields are 
basically given by the formulas for uniform translation: the remain¬ 
ing terms are corrections of order wr/c. Physically, since each star 
moves roughly in a straight path over a time scale shorter than v/a, 
where a = is the acceleration, the fields inside a distance 

cv/a are not affected by time-delay effects, and equal the fields 
due to a dipole in uniform translation. The condition r < cv/a can 
easily be rewritten as 

r < - (28) 

ui 

i.e., it restricts the domain of validity of the approximation to the 
volume inside the light cylinder defined in terms of the orbital an¬ 
gular speed, LJ. Then the electric field generated by star 1 is 

A Bi, (29) 

c 

where Pi is the velocity of star 1. This component of the electric 
field is perpendicular to the magnetic field Bi of the star itself. We 
also add the contribution due to the other star (the one we called 2) 
to the total electric field 

E = E 2 + E 2 ~— ABi + — AB 2 , (30) 

c c 

which, in general, is not perpendicular to the total magnetic field 
B — Bi + B 2 . The component of the induction electric field along 
the magnetic field is 



where p = Pi — P 2 is the stars’ relative speed and h is the unit 
vector along the magnetic field B, in the case of nearly equal-mass 
stars. Recalling that B 2 oc r“^, p oc and the length over 

which E is to be integrated is oc r, we arrive at the approximation 
focr-"/^ 

The meaning of this approximation is the following. Close to the 
point where the two magnetospheres have roughly equal values of 
Bi Ki B 2 ~ B, we find i5|| « VB 2 /C ~ vB/c, i.e. the parallel 
electric field is as big as the maximum achievable. Far from this 
region, i.e. closer to the surface of star 1 (where B 2 <C Bi ~ 
B), the magnetic line will be nearly exactly parallel to Bi, and 
the induction field nearly due to the motion star 1 alone, Ei « 
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vi A Bi/c. thus _E|| « vBijc <C vB/c, i.e. the parallel electric 
field acts as a small perturbation. 

As an order of magnitude, we take for E the value estimated 
above for the in vacuo case: 


7 X 


M 


2.8 X Mq 


1/2 


B2 / Rns \ 


5/2 


1012 0/ \ r J 


eV, 

(32) 


where r = i ?/2 is the half-distance between the stars, B 2 is the 
dipolar component of the magnetic field at the surface of star 2 , 
and 77 « 1 is a fudge factor depending on the relative directions 
of pi and p 2 , the choice of the magnetic field line along which 
we perform the numerical integration, the choice of phase in its 
travel from dayside to nightside and vice This analytical 

computation is borne out by a numerical integration, in fig. for 
arbitrary orientation. The scaling with r applies to all orientations 
we investigated. 

The reason why we chose, among all possible magnetic field 
line, the last one closing onto a star, in the numerical estimate 
above, is that £ has a maximum for this line. In all other cases, 
the total magnetic field is dominated by the nearer star, the induc¬ 
tion electric field is also dominated by the nearer star, but for these 
Ei'Bi = 0. The only effect is due to the component E 2 « VB 2 /c, 
which however decreases due to the larger distance from the star 2 . 
This is shown in fig.[^ where £ is computed for fixed orbital sep¬ 
aration, but for different distances d from star 1 at which the line 
crosses the line joining the two stars’ centers. 


4 CONCLUSIONS 

In this paper, we have begun an investigation of the electrodynam¬ 
ics of binary pulsars. In analogy with IRAPs, we have first deter¬ 
mined the rate of energy loss due to (quadrupolar, not dipolar) ra¬ 
diation, showing that it is dwarfed by GW losses, but may yet be 
detectable due its strong transient character, and its near period¬ 
icity. We have also discussed the presence of transient induction 
electric fields, which, in vacuo, cause extremely large electromo¬ 
tive forces along magnetic field lines in the corotating magneto¬ 
sphere. We have also argued, on the basis of an analogy with good 
conductors, that the skin-depth of these transient fields in a realistic 
(i.e., charge-rich) magnetosphere, is large once the stars are close 
to merging, while it is truly thin at large distances from the star. 
This implies that free charges in the joint magnetosphere will be 
subject to large electric fields with a component parallel to the total 
magnetic field. The consequences of this simple observation will 
be discussed elsewhere. 


and 

(/3 A n)^ = [1 — (cos ujt cos 9 -f sin uut sin 6 cos />)^] . 

(34) 

The time average over the orbital period is now trivial, and we ob¬ 
tain: 

^4/ o 2 rl /•2n 

P = ^JjcosoJ^ .0x 

(/ix cos 9 + fly sin 9 cos cf) + fi^ sin 9 sin 0)^ x 
(1 — ^ cos^ ^ ~ ^ ^ ~ 

= ^(3Mi+4Mz) (35) 
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APPENDIX 

Here we derive eq. [T^ from eq. We take the axis z to be per¬ 
pendicular to the plane of the binary star motion, so that, in Carte¬ 
sian coordinates, P = a;^/3(coscut,sintuf,0). Also, pi — p 2 = 
{fix, fly, fiP), and h = (cos 9, sin 9 cos 4>, sin 9 sin ((>). We thus find 

h ■ {pi — P 2 ) = fix cos 9 + fly sm9cos(l> + fi^. sinS sinc^i , (33) 

^ Strictly speaking, there are also relative orientations of pi and p 2 for 
which the above integral necessarily vanishes. However, this happens only 
when both the magnetic moments are strictly parallel to to, a sufficiently 
rare occurrence that we may neglect in the following. 
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